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Abstract 

We discuss the origin of the difference between the harmonic value of the 
width of the multiphonon giant resonances and the smaller observed value. 
Analytical expressions are derived for both the effective width and the average 
cross-section. The contribution of the Brink- Axel mechanism in resolving the 
discrepancy is pointed out. 
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In a series of publications we have advanced the idea that the excitation of the 

double giant dipole resonances (DGDR) in heavy-ion collisions proceeds in two distinct, in- 
coherent, ways: the usual sequential excitation via the single giant dipole resonance (GDR), 
usually referred to in the literature as the "harmonic" cross section [Q, and the fluctuation 
excitation that involves the "internal" decay of the GDR into the complicated background 
states followed by a collective excitation of a GDR on these states, a manifestation of the 
Brink-Axel mechanism. Of course in a fully microscopic approach to the excitation process, 
the Brink- Axel effects are also contained. However, the usefulness of the approach of [|I],^ is 
the clear identification of the average cross-section with a "direct" excitation and the fluctu- 
ation one with the higher-order "non-direct" contribution. This allows identifying different 
time scales with the excitation process which makes possible an unambiguous treatment of 
the eventual decay of the multiphonon states [2]. The question of the "enhancement" of 
the experimental cross-section with respect to the harmonic value is easily addressed within 
the average+fluctuation model However, the enhancement problem, especially in ^^^Xe 
still remains not fully resolved. There has also been some discussion in recent years related to 
the fact that the observed width of the double giant dipole resonance deviates from the har- 
monic value of twice the single giant resonance width. As we have commented previously, |I| 
this is a natural result of the incoherent contributions to the cross section, even when the 
resonances themselves are purely harmonic. In fact, due to the energy dependence of the 
incoherent contributions to the excitation cross sections, we expect the effective n-phonon 
width to be energy dependent. The purpose of this paper is to derive an analytical expres- 
sion for the effective energy- dependent width of the n-phonon excitation distribution. For 
this purpose we use multistep distorted wave series in conjunction with statistical averaging 
that allows us to significantly extend the result of our previous publication [1] . The method 
enables us to calculate the inclusive n-phonon cross-section in closed form and allows us to 
identify the physical quantity responsible for the cross-section enhancement. This quantity 
is the product of the damping width of the single phonon resonance times the collision time 
evaluated at the grazing impact parameter bo- 
lt is convenient to introduce projection operators: P will stand for the entrance channel, 
di the 1 phonon channel, qi the corresponding fine structure subspace etc. We use multistep 
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distorted wave picture. Thus for the excitation of e.g. 2 phonons, one has 



4-),/ VGY' V 



(+) 



(1) 



In the above V stands for the Coulomb interaction responsible for the excitation, x^"*"' 
stands for the appropriate distorted wave and Gi"*"^ is the Green's operator in the subspace 
di + gi. This is given by 



ydidi ^di "diqi yqiqi 

,(+) 



(2) 



„(+) n^-^' 

where, g'^^^ = {E^+) - - VijG+Vj^)~\ G^^ = = Kei + K and 

E^'^^ — E + ie, with hrei representing the Hamiltonian of relative motion and hi the nu- 
clear Hamiltonian of ion i. The interaction v, couples the different subspaces. 

The interaction V connects the ground state to the one-phonon doorway state di . It can 
also connects a state in gito a one-phonon Brink-Axel state. Thus the amplitude that enters 
into the definition of the cross-section cr^^^ separates into two pieces 
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The fine structure states {qi} are quite complicated many particle-many hole configu- 
rations. In this respect, a statistical treatment of the second, fluctuation, contribution is 
called for. Thus we take the energy average of a^"^^ to represent the theoretical cross-section 
to be compared to the data. When performing the energy average, the cross term in Eq. (4) 



vanishes since it is linear in the assumed random coupling Vq^di ■ Thus a^^^ is given, as usual, 
by the incoherent sum of two terms. The first term is the coherent "harmonic" cross-section. 
Here the Green's function {E — — Vd^q^Gq^Vq^d^Y^ is replaced by an average over the qi 
states [E — Hdi ~ '^diq^Gq^Vq^di) ^- The correction to this approximation involves fiuctu- 
ations that corresponds to the process di — > q'l — > di, which is negligible. The average 
VdiqiGq^Vq^di IS as usual written as A^jj — i-^, which defines the resonance shift A^^ and 
width respectively. These two quantities satisfy a dispersion relation. The fiuctuation 



contribution to cr^^) involves an average over the qi states. 
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This average is easily performed if we recall that the gi Green's function is a large sum 
of random contributions. Thus, writing first 

ki) {(111 



4^^) \Qi) 



/dk 
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where we have employed the scattering eigenstates of hrei, | 
average above into 



(4) 



we can reduce the 



E — £q^— — Ek 



■E 
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The average over {qi} makes the double qi sum collapses into a single sum (gi = q'l)- Fur- 
ther, since there are many qi states, we replace the qi sum by an integral S= / dEg^p (e^J ^ 
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where p{sq^) is the density of the qi states. The Eg integral can be easily performed if we 
assume that the numerator is slowly varying. Only the two Eg — poles will contribute. Thus 
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where we have introduced the spreading width of the doorway state di, defined by 



(6) 



(7) 



and also introduced a representative excited intrinsic state |0 which acts as a ground 

state for the Brink- Axel phonon excitation. 

The energy correlation function |— Fog^ +i(Ek — Sk')\~^ controls the magnitude of crj/^- 
For Ek — and Fq^i — > 0, the contribution is much smaller than af^ since only one k 
integral survive in the former. On the other hand, if Fogj is much larger than the range of 
values of Sk — Sk' which an relevant for the integral, then al/ would be =^ times a regular 

i Oqi 
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cross section. However, Fogi is by definition <^ F^, and therefore this last case does not occur. 
Therefore we may simply say that the correlation function introduces a characteristic time, 
the "correlation" time, which depends on the bombarding energy. Accordingly we write 



(2) 

for cTy/ the following 



4f =r^^re(E) [xy>,l\V\d{s,,)){d^\g,V\Ox, 



(8) 



where the integrals f f 



have been set equal to unity 



since the distorted waves sts {x'^^^} is complete. 

Now since \ 1 (sdj^, contains one Brink- Axel phonons in the final state, where 

as the corresponding amplitude in the harmonic cross-section, \ 2 \V\ 1^, contains 2, we 

(x'-r\i\v\d) , 
can make the approximation \ , . f = -t=. 

{x'f\2\V\l) V2! 



Thus we find for crj^-*, the following reasonable approximate form 



^fi - 2 
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Thus 



(10) 



The second term above gradually becomes insignificant as the bombarding energy in- 
creases. 

The calculation of the cross-section for higher number of phonons follows a similar proce- 
dure. Keeping track of the number of routes that can be followed to reach the final "state" , 
we have for a^^^ 



(11) 



The generalization to n-phonon is straightforward 

^ " ^ n\ n + k\kj 

k=l ^ ' 
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Note that ai"'^ is the cross-section for the excitation of n-phonon states that proceeds 
though n — l,n — 2, phonon states which have finite hfe times (width). In this respect 
we are generaUzing the concept of harmonic cross-section. The consideration of the width 
of the intermediate states and the fluctuations that result from the decay of these states go 
hand in hand. 

Our analytical formula for the average cross-section is quite reasonable when compared 
to numerical solution of the evolution equation for the density matrix of the system reported 
in ref. [3] if the correlation time is identified with the collision time at the grazing impact 
parameter, Tc = K/lv, Fig. 1. In view of this we are confident in applying our theory to 
calculate other observables such as the width of the final channel. 

Taking a Breit-Wigner form for the spectrum of each of the component in (12) and 
evaluating this at the peak, we can define an effective width as 



(n) n-1 (n) /, x 



or 
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^^•^^ = ^-1 . („) (^r,J . (14) 



Eq. (14) for Fg^y (E) shows that the effective width attains the harmonic value nF^^ at 
high energies. At low bombarding energy the value can be significantly smaller than nF^^, 

owing to the Brink- Axel reduction factor I 1 -|- — — ^ . This is demonstrated in 

Fig. 2. The B-A mechanism thus supphes a correlation between the experimental effective 

width and the corresponding cross-section. This is an important result which further helps 

(2) i+^^sXI 

shed light on the excitation mechanism. For two phonons we have ^ = -p^. This 
relation is independent of the system. In figure (3) we exhibit the available data taken 

(2) (2) 

from ref. [6]. In presenting the data, we have changed the variable into with 

'^harm being the value of (Tc when the width of the intermediate one-phonon state is set 
equal to zero. Although it is claimed that F^^^^ = \/2Fi, our result does not exclude a more 
subtle connection between F^^^ and Fi. 
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The above discussion suggests defining an effective phonon 



number ^/Z i+"j:'^kjr,) which is smaller than n. This gives an immediate quali- 
^ 

tative explanation for the enhancement in the cross-section too. Since rig// < n, the energy 
^Ueff < nEi. Accordingly, if one were to use damped anharmonic oscillator model with 
the following sequence of energies, — > — > 2Ei — j^^Ei — > 3Ei — -^^^^^^^^Ei, one 
would obtain larger cross-sections cr^^^; a^^\ ... when compared to ai^\ ai^\... In a way, this 
supplies a simple connection between the direct + fluctuation model and the anharmonic 
models [7-9]. In Fig. 4 we show the result of a calculation following this line. The multi- 
phonon cross-section calculated for the effective damped anharmonic oscillator (long-dashed 
line), whose spectrum is given by -Ene// — ^e// -E'l- The semiclassical coupled channels 
model [4b] is used to generate this result. The widths of the intermediate states are just 
the "harmonic" ones Tn = nTi. Also shown is the cross-section from Ref. [3] (solid line) 
and the simple undamped harmonic oscillator cross-section also calculated according to [4b] 
(small dashed fine). The results shown are for the system ^°^Pb-|-^°^Pb. The agreement 
of our damped anharmonic oscillator, coupled channels calculation with those of Ref. [3], 
which is based on the solution of an evolution equation for the excitation probability with 
loss and gain terms, is excellent. 

In conclusion, the direct-|-fluctuation model (Brink- Axel mechanism) supplies a natural 
framework to discuss the cross-section enhancement and the width reduction of the multi- 



phonon states. The enhancement factor of the cross-section and the width reduction 

p(n) 



can be co-related by the effective number of phonons. This concept should be quite use- 



ful in future work on multiphonon physics. Further, our results here should also be relevant 
to other multistep reaction phenomena such as preequilibrium processes [10]. 
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FIGURE CAPTIONS 



Fig. 1: The quantity as a function of the bombarding energy for the ^^^Ph + 208p|-j 
system. 

Fig. 2: The effective observed width r^^j, Eq. (14) vs. the bombarding energy for the 
208pb+208pb system. Clearly T^^j^ = T (di), and it is set equal to 4.0MeV. See text for 
details. 

Fig. 3: The effective width reduction vs. the cross-section enhancement. The data are 

(2) 

from Ref. [6]. The "data" for the cross-section ratios is which is just the values 

of Ref. [6], ^f^, multiplied by ^\tf^\ since 4aL. = ^^c^ (^di = 0). See text for 
details. 

Fig. 4: The multiphonon cross-section calculated for the effective damped anharmonic os- 
cillator (long-dashed line), whose spectrum is given by En^^^ — ng// E^. The semi- 
classical coupled channels model [4b] is used to generate this result. The widths of 
the intermediate states are just the "harmonic" ones r„ = nVi. Also shown is the 
cross-section from Ref. [3] (solid hne) and the simple undamped harmonic oscillator 
cross-section also calculated according to [4b] (small dashed line). The results shown 
are for the system ^°^Pb-|-^°*Pb. See text for details. 
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